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We show that the special relativistic dynamics when combined with quantum mechanics and 
the concept of superstatistics can be interpreted as arising from two interlocked non-relativistic 
stochastic processes that operate at different energy scales. This interpretation leads to Feynman 
amplitudes that are in the Euclidean regime identical to transition probability of a Brownian particle 
propagating through a granular space. Some kind of spacetime granularity could be held responsible 
for the emergence at larger scales of various symmetries. For illustration we consider also the 
dynamics and the propagator of a spinless relativistic particle. Implications for doubly special 
relativity, quantum field theory, quantum gravity and cosmology are discussed. 

PACS numbers: 03.65.Ca, 03.30.+p, 05.40.-a, 04.60.-m 
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Introduction. — The concept of "emergence" plays an 
important role in quantum field theory and, in partic- 
ular, in particle and condensed matter physics, since 
it embodies the essential feature of systems with sev- 
eral interlocked time scales. In these systems, the ob- 
served macroscopic-scale dynamics and related degrees 
of freedom differ drastically from the actual underlying 
microscopic-scale physics Superstatistics provides a 
specific realization of this paradigm: It predicts that the 
emergent behavior can be often regarded as a superposi- 
tion of several statistical systems that operate at differ- 
ent spatio-temporal scales ||. In particular, many ap- 
plications have recently been reported, in hydrodynamic 
turbulence 0] , turbulence in quantum liquids [H , pattern 
forming systems Q or scattering processes in high-energy 
physics [7[. 

The essential assumption of the superstatistics scenario 
is the existence of sufficient spatio-temporal scale separa- 
tions between relevant dynamics within the studied sys- 
tem so that the system has enough time to relax to a lo- 
cal equilibrium state and to stay within it for some time. 
In practical applications one is typically concerned with 
two scales. Following Ref. Q, we consider an intensive 
parameter £ that fluctuates on a much larger time scale 
than the typical relaxation time of the local dynamics. 
The random variable £ can be in practice identified, e.g., 
with the inverse temperature 0, Q , friction constant || , 
volatility Q or einbein (Io| . On intuitive ground, one 
may understand the superstatistics by using the adia- 
batic Ansatz. Namely, the system under consideration, 
during its evolution, travels within its state space X (de- 
scribed by state variable A <E X) which is partitioned into 
small cells characterized by a sharp value of (. Within 
each cell, the system is described by the conditional dis- 
tribution p{A\C,). As ( varies adiabatically from cell to 
cell, the joint distribution of finding the system with a 



sharp value of C in the state A is p(A,Q = p(A\C)p(() 
(Baycs theorem). The resulting macro-scale (emergent) 
statistics p(A) for finding system in the state A is ob- 
tain by eliminating the nuisance parameter £ through 
marginalization, that is 



P (A) 



p(A\C)p(C) d(. 



(1) 



Interestingly enough, the sufficient time scale separation 
between two relevant dynamics in a studied system al- 
lows to qualify superstatistics as a form of slow modula- 
tion [m. 

In this Letter, we recast the Feynman transition am- 
plitude of a relativistic scalar particle into a form, which 
(after being analytically continued to imaginary times) 
coincides with a superstatistics marginal probability (flj . 
The derivation is based on the Levy-Khinchine theo- 
rem for infinitely divisible distributions (l2l . 13 [, and 



for illustration we consider the dynamics and the prop- 
agator of a Klein-Gordon (i.e., neutral spin— 0) parti- 
cle. Our reasonings can be also extended to charged 
spin— 0, spin— i, Proca's spin— 1 particles and to higher- 
spin particles phrased via the Bargmann— Wigner wave 
equation (lol | . Further generalization to external electro- 
magnetic potential has been reported in Refs. [H[i3. 

We also argue that the above formulation can be looked 
at as if the particle would randomly propagate (in the 
sense of Brownian motion) through an inhomogeneous 
or granular medium ( "vacuum" ) [l4| . 

Our argument is based upon a recent observation [^, 
Hot 14 1 that the Euclidean path integral (PI) for rel- 



ativistic particles may be interpreted as describing a 
doubly-stochastic process that operates at two separate 
spatio-temporal scales. The short spatial scale, which is 
much smaller than particle's Compton length A c = l/mc 
(h= 1), describes a Wiener (i.e., Galilean relativity) pro- 
cess with a sharp (Galilean-invariant) Newtonian mass. 
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The large spatial scale, which is of order A c , corresponds 
to distances over which the fluctuating Newtonian mass 
changes appreciably. At scales much larger than \ c the 
particle evolves according to a genuine relativistic dy- 
namics, with a sharp value of the mass coinciding with 
the Einstein rest mass. Particularly striking is the fact 
that when we average the particle's velocity over the 
structural correlation distance (i.e., over particle's A c ) 
we obtain the velocity of light c. So the picture that 
emerges from this analysis is that the particle (with a 
non-zero mass!) propagates over the correlation distance 
A c with an average velocity c, while at larger distance 
scales (i.e., when a more coarse grained view is taken) 
the particle propagates as a relativistic particle with a 
sharp mass and an average velocity that is subluminal. 
Quite remarkably, one can observe an identical behavior 
in the well-known Feynman's checkerboard PI III 17 1 to 
which the transition amplitude ([I]) reduces in the case of 
a relativistic Dirac fermion in 1 + 1 dimensions T3, 3 ■ 
A considerably expanded presentation including the 
issue of reparametrization invariance, bibliography, and 
proofs of the main statements and formulas is given in a 
companion paper 



p(x',t'\x,t) satisfies the ESE ©, was found in Ref. [§]. 
The key is to note that in order to have ([2]) satisfied by p, 
the rescaled PDF w((, T) = w{C/T, T)/T should satisfy 
the ESE for homogeneous Markovian process 

w(C,h+t 2 ) = [ dC'«,(C',*iMC-C',*2). (4) 

Jo 

Consequently the Laplace image fulfills the functional 
equation with *i,t2 £ R + - By assuming continuity in T, 
it follows that the multiplicative semigroup u>(p c , T)t>o 
satisfies w(p ( ,T) = {w(p., 1)} T - From this we see that 
the distribution of £ at T is completely determined by 
the distribution of f at T = 1. In addition, because 
w(p ( ,l) = {w(p c ,l/n)} n for any n <E N + , w(C, 1) is in- 
finitely divisible. The Levy-Khinchine theorem 



then ensures that log w(p ( . 
generic form 



T) = -TF(p ) must have the 



logw(p ( ,T) = -T [ a P( + ^ (1 - e- p < u )v(du) ) , (5) 



14] 



where a > is a drift constant and v is 
some non-negative measure on (0, oo) satisfying 
/ R+ min(l, u)v{du) < oo. Finally the Laplace inverse 
of w(p ( ,T) yields uj((,T). Once uj(( 7 T) is found, then 
p(x',t'\x,t) possesses a PI representation on its own. 
What is the form of the new Hamiltonian? To this end 
dzp(y,t"\z,t')p(z,t'\x,t) , (2) we rewrite (J3]) in Dirac operator form as @ 



Superstatistics path integrals. — When a conditional 
probability density function (PDF) is formulated through 
a PI, then it satisfies the Einstein- Smoluchowski equation 
(ESE) for continuous Markovian processes, namely fl5j ] 



p(y,t"\x,t) 



with t' being any time between t" and t. Conversely, 
any transition probability satisfying ESE possesses a PI 
representation [16[ . In physics one often encounters prob- 
abilities formulated as a superposition of Pi's, 

p(x',t'\x,t) 

dCw(C,T) f ( ) [dxdp] e J? dT(ipx-CH(p,x)) 
Jx(t)=x 

(3) 

Here u>(C,T) with T = t' - t is a normalized PDF de- 
fined on K + x R + . The form typically appears in 
non-perturbativc approximations to statistical partition 
functions, in polymer physics, in financial markets, in 
systems with reparametrization invariance, etc. The ran- 
dom variable Q is then related to the inverse temperature, 
coupling constant, volatility, vielbein, etc. 

The existence of different time scales and the flow of 
the information from slow to fast degrees of freedom typi- 
cally creates the irreversibility on the macroscopical level 
of the description. The corresponding information thus 
is not lost, but passes in a form incompatible with the 
Markovian description. The most general class of distri- 
butions uj((, T) on R + x R + for which the superposition of 
Markovian processes remain Markovian, i.e., when also 



p(x',t'\x,t) = (x'\ / dtw(t,T)e-< H \x) 



-TF(H)\ 



{x'\{w(H,l)} T \x) 



(0) 



Hence, the identification H(p, x) = F(H(p, x)) can be 
made. Here one might worry about the operator-ordering 
problem. For our purpose it suffices to note that when H 
is ^-independent, the former relation is exact. In more 
general instances the Weyl ordering is a natural choice 
because in this case the required mid-point rule follows 
automatically and one does not need to invoke the gauge 
invariance [9|, l29( . In situations when other non-trivial 
configuration space symmetries (such as non-holonomic 
symmetry) are required, other orderings might be more 
physical [9(. 

Emergent Special Relativity. — The Feynman transi- 
tion amplitudes (or better its Euclidean version — tran- 
sition probabilities) naturally fits into the structure of 
superstatistics Pi's discussed above. 

Note first that the choice a = and i/(du) = 



l/(2^u 3 ' 2 )du leads to F(p ( ) = y /p^. This identi- 
fies w(£, T) with the (unshiftcd) Levy distribution with 
the scale parameter T 2 /2. Moreover, when H(p, x) = 
p 2 c 2 + m 2 c 4 then ([3]) can be cast into the form (see also 
Refs. jiESIll) 
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n X {T)=x' ( nT 

I [dx dp] cxp < / dr 

'x(0) = x \Jo 



ip ■ x — cy/p 2 + m 2 c' 



dm fi (m,Tc 2 ,Tc 



2 v^ m *\ 



x{T) = x' ( pT 

[dx dp] cxp < / dr 
x(0) = x [Jo 



. . p' 

ip ■ x — — — mc 
v 2m 



r 



(7) 



where t' — t = T, and 

{a/b)P/ 2 



i p (z,a, b) 



,P-1 p -(az+6/^)/2 



2K p (Vab) 



(8) 



is the generalized inverse Gaussian distribution [13j (.ffp 
is the modified Bessel function of the second kind with 
index p). The LHS of ([7]) represents the PI for the free 
spinless relativistic particle in the Newton- Wigner repre- 
sentation [lH . The full Klein-Gordon (KG) kernel which 
also contains the negative-energy spectrum can be ob- 
tained from ([7]) by considering the Fcshbach-Villars rep- 
resentation of the KG equation and making the substi- 
tution El 



fi (m,tc 2 ,tc 2 m 2 ) 



1 



i — ^ ■ 



-^^h(m,\t]c 2 ,]tM.(9) 



The matrix nature of the smearing distribution (173 
is the Pauli matrix) naturally includes the Feynman- 
Stuckelberg causal boundary condition and thus treats 
both particles and antiparticles in a symmetric way 
fioj ] . When the partition function is going to be calcu- 
lated, the trace will get rid of the sgn(i) term and 1/2 is 
turned to 1. 

The explicit form of the identity ([7]) indicates that m 
can be interpreted as a Galilean-invariant Newton-like 
mass which takes on continuous values distributed ac- 
cording to f 1 (m, Tc 2 ,Tc 2 m 2 ) with (to) = rn+l/Tc 2 and 
var(m) = m/Tc 2 + 2/T 2 c 4 . Fluctuations of the Newto- 
nian mass can be then depicted as originating from parti- 
cle's evolution in an inhomogeneous or granular medium. 
Granularity, as known, for example, from solid-state sys- 
tems, typically leads to corrections in the local disper- 
sion relation [20( and hence to alterations in the local 
effective mass. The following picture thus emerges: on 
the short-distance scale, a non-relativistic particle can 
be envisaged as propagating via classical Brownian mo- 
tion through a single grain with a local mass m. This 
fast-time process has a time scale ~ 1/mc 2 . An aver- 
aged value of the local time scale represents a transient 
temporal scale (1/mc 2 ) = 1/mc 2 which coincides with 
particle's Compton time T c — i.e., the time for light to 
cross the particle's Compton wavelength. At time scales 
much longer than T c (large-distance scale), the proba- 
bility that the particle encounters a grain which endows 
it with a mass m is f 1 (m, Tc 2 , Tc 2 m 2 ) . As a result one 



may view a single-particle relativistic theory as a single- 
particle non-relativistic theory where the particle's New- 
tonian mass m represents a fluctuating parameter which 
approaches on average the Einstein rest mass to in the 
large t limit. We stress that t should be understood as the 
observation time, a time after which the observation (po- 
sition measurement) is made. In particular, during the 
period t the system remains unperturbed. One can thus 
justly expect that in the long run all mass fluctuations 
will be washed out and only a sharp time-independent 
effective mass will be perceived. The form of (m) identi- 
fies the time scale at which this happens with t ~ 1/mc 2 , 
i.e. with the Compton time T a , It should be stressed 
that above mass fluctuations have nothing to do with the 
Zitterbewegung which is caused by interference between 
positive- and negative-energy wave components. In our 
formulation both regimes are decoupled. 

We may also observe that by coarse-graining the ve- 
locity over the correlation time T c we have 



(N) 



m 



(10) 



So on a short-time scale of order A c the spinless relativis- 
tic particle propagates with an averaged velocity which 
is the speed of light c. But if one checks the particle's 
position at widely separated intervals (much larger than 
A c ), then many directional reversals along a typical PI 
trajectory will take place, and the particle's net velocity 
will be then less than c — as it should be for a massive 
particle. The particle then acquires a sharp mass equal 
to Einstein's mass, and the process (not being hindered 
by fluctuating masses) is purely Brownian. This conclu- 
sion is in line with the well-known Feynman checkerboard 
picture [14|,|l_7| to which it reduces in the case of (1 + 1)D 
relativistic Dirac particle. 

Robustness of emergent special relativity. — Under- 
standing the robustness of the emergent Special Rela- 
tivity under small variations in the mass-smearing distri- 
bution function f 1 can guide the study of the relation be- 
tween Einsteinian SR and other deformed variants of SR, 
such as Magueijo-Smolin and Amelino-Camelia's doubly 
special relativity [2ll . |22j , or (quantum) K-Poincare defor- 
mation of relativistic kinematics (23)]. In DRS models a 
further invariant scale I is introduced, besides the usual 
speed of light c, and I is typically considered to be of 
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order of the Planck length. A small variation Si i of the 

2 

smearing function originates the new Hamiltonian 



H = \ + (l + |) ^P 2 c 2 + mV + | , (11) 

with ei = -2(1 + e /2)^ej (sec Ref. [3 for details). 
By setting 



£i=2 



1 



1 — c 2 m 2 £ 2 



1 , 6 2 



4c 6 m 4 ^ 2 
1 — c 2 m 2 £ 2 



the new Hamiltonian 77 can be easily identified with 



2 c 2 ^ =p \/p 2 (l — m 2 c 2 £ 2 ) + m 2 c 2 
1 — m 2 c 2 £ 2 



(12) 



which coincides with the Magueijo-Smolin's doubly spe- 
cial relativistic Hamiltonian, in, say, its version [24j . It 
should be stressed that the Hamiltonian (|11|) (when also 
negative energy states are included) violates CPT sym- 
metry. This is a typical byproduct of the Lorentz sym- 
metry violation in many deformed SR systems. 

For the Hamiltonian (TT2|) a relation analog to ([7]) 
holds, where now the smearing function has the form 
f i (m, Tc 2 A, Tc 2 m 2 \) with A = 1/(1 - m 2 c 2 £ 2 ). The cor- 
relation distance 1 / mcA can be naturally assumed as the 
minimal size £ GRAIN of the "grain of space" of the poly- 
crystalline medium, which is linked to the new invariant 
scale £ by 



1 



mcA 



= A c (l -m 2 c 2 £ 2 ). 



(13) 



By tuning the size £ GRAIN of these "grains of space" it 
is possible to pass continuously from Lorenz symmetry 
to other different symmetries, as those enjoyed by DSR 
models. We can in principle speculate that each large 
scale symmetry could originate from a specific kind of 
space (time) granularity. 

Quantum field theory. — The superstatistics transition 
probability © was constructed on the premise that H is 
associated with a single particle. Of course, a single- 
particle relativistic quantum theory is logically unten- 
able, since a multi-particle production is allowed when- 
ever the particle reaches the threshold energy for pair 
production. In addition, Leutwyler's no-interaction the- 
orem (25| prohibits interaction for any finite number of 
particles in the context of relativistic mechanics. To 
evade the no-interaction theorem it is necessary to have 
an infinite number of degrees of freedom to describe inter- 
action. The latter is typically achieved via local quantum 
field theories (QFT). 

It should be underlined in this context that the PI for 
a single relativistic particle is still a perfectly legitimate 
building block even in QFT. Recall that in the standard 
perturbative treatment of, say, generating functional for 



a scalar field each Feynman diagram is composed of in- 
tegrals over product of free correlation functions (Feyn- 
man's correlators): 

1 [°° 

A F (y 1 ct y ;z,ct z ) = - / drsgn(r - ty)p(y,r\z,t z ), (14) 
^ J — 00 

and may thus be considered as a functional of the PI 
p(x', i'|x, t). In fact, QFT in general, can be viewed as 
a grand-canonical ensemble of fluctuating particle his- 
tories (worldlines) where Feynman diagrammatic repre- 
sentation of quantum fields depicts directly the pictures 
of the world-lines in a grand-canonical ensemble. This is 
the so-called "worldline quantization" of particle physics, 
and is epitomized, e.g., in Feynman's worldline represen- 
tation of the one-loop affective action in quantum elec- 
trodynamics (26[, in Kleinert's disorder field theory [2^ 
or in the Bcrn-Kosower and Strassler "string-inspired" 
approaches to QFT [28^ 

Because of HH), the relationship between bosonic 
Bern-Kosower Green's function Gb(ti,T2) and the PI 
p(x', i'|x, t) can be found easily through the known func- 
tional relation between Gb and Ap, cf. Refs. (28j |. 

Gravity and Cosmology. — When spacetime is curved, 
a metric tensor enters in both Pi's in (7) in a different 
way, yielding different "counterterms" fisl [29j . For in- 
stance, in Bastianelli-van Nieuwenhuizen's time slicing 
regularization scheme [29I ] one has (when h is reintro- 
duced) 



t_ 9 ij PiPj ^ „ ; 

2m * 2m + 8m ( ^ + 5 1 « 1 W> 



. r~ h 2 

Vp 2 + m 2 c 2 1 y ygvpiPj + —{R + g Ij r™rj„j + m 2 c 2 

+ h 4 <f>(R,dR,d 2 R) + 0(h e ), (15) 

where g*i , R, T ] kl and $(...) arc the (space-like) pull- 
back metric tensor, the scalar curvature, the Christoffcl 
symbol, and non- vanishing function of R and its first and 
second derivatives, respectively. This causes the super- 
statistics identity (7) to break down, as can be explic- 
itly checked to the lowest order in h. The respective 
two cases will thus lead to different physics. Because 
the Einstein equivalence principle requires that the local 
spacetime structure can be identified with the Minkowski 
spacetime possessing Lorentz symmetry, one might as- 
sume the validity of (0 at least locally. However, in 
different space-time points one has, in general, a differ- 
ent typical length scale of the local inertial frames, de- 
pending on the gravitational field. The characteristic size 
of the local inertial (i.e. Minkowski) frame is of order 
1/lifl 1 / 4 where K = R a 0~ ( sR a ^ 5 is the Kretschmann 
invariant and R a p-ys is the Riemann curvature. Rela- 
tion ([7]) tells us that the special relativistic description 
breaks down in regions of size smaller than A c . For 
curvatures large enough, namely for strong gravitational 
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fields, the size of the local inertial frame can become 
smaller than A c , that is 1/|^C| 1//4 < A c . In such regions 
the special relativistic description is no more valid, and 
according to (O must be replaced by a Newtonian de- 
scription of the events. For instance, in Schwarzschild 
geometry we have K = 12 r 2 /r 6 , and the breakdown 



should be expected at radial distances 



< 



1/3 



(r s is the Schwarzschild radius) which are — apart from 
the hypothetical case of micro-black holes (where A c ~ 
r s ) — always deeply buried below the Schwarzschild 
event horizon. In the cosmologically relevant Fricdmann- 
Lcmaitre-Robertson- Walker (FLRW) geometry, we have 
K = 12 (d 4 + a 2 a 2 ) / (ac) 4 , and the breakdown happens 
when (d 4 + a 2 d 2 ) > (ac/A c ) 4 , where a{t) is the FLRW 
scale factor of the Universe and d = da /At. Applying 
the well-known Vilenkin-Ford model [3(j for inflationary 
cosmology, where a(t) is given by: a(t) = A^smh(Bt) 
with B = 2cy/ A/3 (A is the cosmological constant), we 
obtain a temporal bound on the validity of local Lorentz 
invariance, which, expressed in FLRW time, is 



t < — arcsinh 
~ B 



BK 



(8c 4 -(BX c ) 4 ) 



4\V 4 



t. (16) 



By using the presently known |3l( value A ~ 10~ 52 m -2 
and the r-lepton Compton's wavelength \ T c ~ 6.7 x 
10~ 16 m (yielding the tightest upper bound on t), we ob- 
4 x lCT 24 s. Note that, since B\ c <C c, then 



A c /c = t a . Such a violation of the local Lorentz in- 



tain t 
i 

variance naturally breaks the particle-antiparticle sym- 
metry since there is no unified theory of particles and 
antiparticles in the non-relativistic physics — formally 
one has two separate theories. If the resulting matter- 
antimatter asymmetry provides a large enough CP asym- 
metry then this might have essential consequences in the 
early Universe, e.g., for leptogenesis. In this respect, t is 
compatible with the nonthermal leptogenesis period that 
typically dates between 10 _26 -10 _12 s after the Big Bang. 

Conclusions and perspectives. — The new superstatis- 
tics PI representation of a relativistic point particle in- 
troduced in this Letter, realizes an explicit quantum me- 
chanical duality between Einsteinian and Galilean rela- 
tivity. It also makes explicit how the SR invariance is 
encoded in the grain smearing distribution. Notably, the 
exact LS of a spacetime has no fundamental significance 
in our analysis, as it is only an accidental symmetry of 
the coarse-grained configuration space in which a par- 
ticle executes a standard Wiener process. In passage 
from grain to grain particles's Newtonian mass fluctu- 
ates according to an inverse Gaussian distribution. The 
observed inertial mass of the particle is thus not a funda- 
mental constant, but it reflects the particle's interaction 
with the granular vacuum (cosmic field). This, in a sense, 
supports Mach's view of the phenomenon of inertia. 

Interactions can be included in our framework in two 
different ways. The interaction with a background field 



(such as electromagnetic field) can be directly treated 
with the superstatistics prescription (|7j) , see [10j . On the 
other hand, the multi-particle interactions can be con- 
sistently formulated by "embedding" the relativistic PI 
in QFT via the worldline quantization. Such an em- 
bedding may help to study several cosmological impli- 
cations of systems with granular space. If any of such 
systems quickly flows to the infrared fixed point, any di- 
rect effect due to the space discreteness, and related SR 
violation, might be insignificant on cosmological scales 
(where Lorentz and diffcomorphism invariance are re- 
stored), while it might be crucial in the early Universe, 
e.g., for leptogenesis and the ensuing baryogenesis. Con- 
sequences on the detailed structure of the Cosmic Mi- 
crowave Background spectrum will be explored in future 
work. 

The presented approach implies a preferred frame. In 
this connection it is worth of noting that, despite the fact 
that ([7]) is not manifestly LS invariant, one may use the 
Stiickelberg trick and introduce a new fictitious variable 
into the PI (0, in such a way that the new action will 
have the reparametrization symmetry, but will still be 
dynamically equivalent to the original action. For rele- 
vant details see Ref. [l4j]. By not knowing the source, one 
may then view this artificial gauge invariance as being a 
fundamental or even a defining property of the relativistic 
theory. One might, however, equally well, proclaim the 
"polycrystalline" picture as being a basic (or primitive) 
edifice of SR and view the reparametrization symmetry 
as a mere artefact of an artificial redundancy that is al- 
lowed in our description. It is this second view that we 
favored here. 

The presented scenario cannot directly accommodate 
the massless particles such as photons (identity (JT)) holds 
true only for m =/= 0). One possibility would be to use the 
PI representation of Polyakov- Wheeler for massless par- 
ticles and try to construct a similar superstatistics dual- 
ity as in the case of massive particles. This procedure is, 
however, not without technical difficulties and currently 
is under investigation. Conceptually is far more simpler 
to assume that the photon has a small mass. At present, 
there are a number of experimental limits to the mass of 
the photons [32|. For instance, tests based on Coulomb's 
law and the galactic vector potential set the upper limit 
of to 7 < 10~ 18 cV/c 2 ~ 10~ 57 g. This gives the domain 



correlation distance for the photon ~ 1/m^ 



10 



43, 



which is bigger than the radius of observable Universe 
(~ 10 26 m) and so in this picture the photon mass does 
not fluctuate — it is a quasi-invariant. 

Finally, this approach should reinforce the links be- 
tween superstatistics paradigm and the approach to 
quantum gravity based on stochastic quantization (33| . 
In particular, the outlined granular space could be a nat- 
ural model for the noise terms in a Parisi-Wu stochastic- 
like quantization approach to gravity. 
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